Phonons in an Inhomogeneous Continuum: Vibrations of an Embedded Nanoparticle 
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The spectrum of acoustic vibrational modes of an inhomogeneous elastic continuum are analyzed 
with application to a spherical nanoparticle embedded in an infinite glass block. The relationship 
of these modes to the discrete vibrational spectrum of a free sphere is studied. The vibrational 
modes of a sphere with a fixed surface are relevant in some situations. Comparisons are also made 
to calculations of mode frequency and damping based on complex valued frequency. 
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I. INTRODUCTION 

In 1882, Lambi calculated the frequencies of vi- 
bration of a free isotropic homogeneous continuous 
sphere. Twenty years agOjS* 3 - this model was applied 
to nearly spherical atomic clusters a few nanometers 
in size (nanoparticles), whose vibrations can be ob- 
served by low frequency inelastic light scattering (Raman 
or Brillouin^ setups depending on the actual particle 
size) and also by femtosecond pump-probe experiments. 6 
However, with exceptions jSiSuLiii the free sphere model 
hardly seems applicable because the samples investigated 
are doped glasses, i.e. nanoparticles solidly embedded 
within a glass matrix. 

Variations on Lamb's calculation attempt to take into 
account the mechanical properties of the matrix. The 
nanoparticle surface can be rigidly fixed^i Assuming out- 
going travelling waves in the matrix, complex valued 
mode frequencies can be found>i2ii£ These approaches 
give discrete sets of vibrational modes, unlike a macro- 
scopic object in which there is a continuum. 

Conversely, the nanoparticle-matrix interface influence 
on vibrations can be ignored^ and bulk phonon modes 
can be used as if the nanoparticle surface were acous- 
tically invisible. This latter approach can be workable 
when other mechanisms, such as size quantization of 
nanoparticle electron energy levels (excitons), provide 
the dominant contribution to the qualitative features of 
light scattering spectra. 

None of the foregoing provides a satisfactory descrip- 
tion of a nanoparticle-matrix system even if they can 
provide valuable insight. What is needed is a calcula- 
tion that properly models the effect of the presence of 
the nanoparticle on the phonon modes. In this paper, we 
present such an approach based on recent calculations on 
core-shell systems M> 

Apart from some general observations made in section 
ITT1 the nanoparticle is idealized as an isolated spherical 
homogeneous isotropic classical continuum object em- 
bedded in similarly idealized glass. Departures from this 
won't be addressed in this paper even though they can be 



experimentally important. Taking into account nanopar- 
ticle size distribution is straightforward. Anisotropy has 
been discussed in a previous paper, 13 and here we use di- 
rectionally averaged sound speeds calculated as described 
therein. Non-sphericity is a more complex problem but 
arbitrary shapes can be investigated using molecular dy- 
namics simulation^ Application of continuum elastic 
boundary conditions at the interface between nanopar- 
ticle and matrix is also an idealization that needs to be 
justified with attention to details at the microscopic level. 
For example, such conditions are generally incorrect for 
epitaxial interfaces^ 

Mode frequencies predicted for a free sphere corre- 
spond to dampened motions at the same frequency of the 
nanoparticle in a sufficiently soft and light matrix. Some 
qualitatively new modes appear corresponding to semi- 
rigid libration (angular oscillation) and rattling (transla- 
tional oscillation) of the nanoparticle as a result of the 
elastic restoring force provided by the matrix. Mode fre- 
quencies for a sphere with a fixed boundary correspond 
to the dampened motions of a nanoparticle in a harder 
and denser matrix. 



II. GENERAL INHOMOGENEOUS SYSTEM 

Although subsequent sections of the paper deal with 
isotropic and homogeneous materials, it is useful to first 
consider a continuous elastic medium which may be in- 
homogeneous and anisotropic. v(R, t) denotes real space 
coordinates of the point of matter which has material 
coordinate R. Deviation from equilibrium is described 
in terms of the displacement field u(R, t) = r(R, t) — R. 
The mass density of the medium is p(R). Note that this 
density refers to material coordinates. The density in 
real space coordinates varies as acoustic vibrations in the 
medium cause local expansion and contraction. 

The elastic constants of the material are Cijki{R). 
Stress (Tij is related to the displacement field through 
o~ij = CijkiUk,i where a summation is implied over re- 
peated indices and the comma implies a partial deriva- 
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tive. The equation of motion is p(R)iii — (Jij.j. 

Six different models, with abbreviations listed in Ta- 
ble [I] have been studied. For some of them (BSM,FSM 
and CSM) u is real valued. For the others, u is complex 
valued, but the imaginary part of u is to be ignored. 

The equation of motion has normal modes u n (R) with 
frequency u„ where n is a mode index. The angular 
frequency to (in rad/s) is related to the wavenumber v 
(in cm -1 ) through v = uj/(200ttc) where c is the speed 
of light (in m/s). Modes are orthonormalized using 



p{R)ui{R) ■ u 3 {R)d 3 R = b lj M p + m 



(1) 



for any two modes i and j, where M p+m is the mass of 
the system (nanoparticle and matrix) . The integration is 
over the volume of the system. Equation JJJ is justified 
in Appendix 1X1 A general disturbance of the system has 
displacement field 

u(R, t) = y](x n sin(u n t) + y n cos(uj n t))u n (R) (2) 



On the basis of the normalization chosen, u n (R) is di- 
mensionless. Therefore x n and y n have units of metres. 

Inserting J5J into (|A1|) and (|A2|) and simplifying using 
(|A6|) and l|A7(l . the total energy of the system is 
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E = ^ M P+m U U X l + Vn 



(3) 
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III. SPHERICALLY SYMMETRIC MODELS 



FSM Free Sphere Model 

BSM Bound Sphere Model 

CFM Complex Frequency Model 

RCFM Rigid CFM 

VCFM Void CFM 

CSM Core-Shell Model 

TABLE I: Abbreviations used throughout this paper. 



of such a mode (a plot of scattered light intensity versus 
wavenumber v) a broadened peak is observed with its 
center at v = Re(w)/(2007rc) and with a half width at 
half maximum (HWHM) of Im(cj)/(2007rc). 

Other worker o 18 ! 19 used the "absolute value of the 
spherical Hankel function" for the outgoing wave, artifi- 
cially reducing the mathematical problem to real-valued 
frequencies. 

There are two important limiting cases of CFM: In the 
rigid complex frequency method (RCFM) the nanoparti- 
cle is made infinitely rigid and heavy. In the void complex 
frequency method (VCFM) the nanoparticle is replaced 
by an empty void. 

CFM, RCFM and VCFM necessarily have complex- 
valued u which blows up exponentially with R and can- 
not be orthonormalized. Equations QJ to 0| and the 
discussion in Appendix |X| are for FSM, BSM and CSM 
only and do not apply to CFM, RCFM and VCFM. 



In this and later sections the displacement field de- 
pends on spherical material coordinates R = (R,9,(ft). 
The nanoparticle ("p") is idealized as a homogeneous and 
isotropic sphere of radius R p and density p p with speeds 
of sound v p L and v?p. 

In the free sphere modeli (FSM) , the nanoparticle sur- 
face at R = R p is free. Mode frequencies are ^f^q 
where q can be spheroidal (SPH) or torsional (TOR), 
I and m are angular momentum and its z-component, 
and n is the harmonic index. For convenience, these 
modes will be called "pseudomodes" when referring to 
situations of a sphere that is only weakly coupled to its 
surroundings. Several incorrect calculations of the FSM 
(SPH,Z=0) modes appear in the literature^ 

In the bound sphere modelii (BSM), the nanoparticle 
surface is rigidly fixed. Mode frequencies are ^^n^ ■ 

In the complex frequency modeliSii^ (CFM), the 
nanoparticle is surrounded by a homogeneous and 
isotropic matrix (" m" ) of density p m and speeds of sound 
v™ and w™. The matrix extends to infinity. Continuity 
of u and force-balance apply at the nanoparticle-matrix 
interface. The boundary condition at large R is that 
u is an outgoing travelling wave. The mode frequencies 
^imnq are complex valued. In spectroscopic observations 



IV. CORE-SHELL MODEL 

Adapting the core-shell modeli^ (CSM) we consider 
a nanoparticle as in the previous section concentrically 
within a much larger homogeneous and isotropic sphere 
(representing the glass matrix) with radius R m . In the 
macroscopic limit R m goes to infinity, and the frequency 
of every vibrational mode varies as 1 / R m and therefore 
goes to zero. Going from the FSM to an embedded sphere 
with a macroscopic matrix means changing from a dis- 
crete set of vibrational modes to a continuum. 

We are primarily interested in Raman or Brillouin ex- 
periments in which light is scattered from the nanopar- 
ticle rather than from the transparent matrix material. 
Low frequency vibrations of the nanoparticle allow the 
scattered photon to have wavenumber a few cm -1 above 
or below those in the incident laser beam. The general 
idea behind this core-shell approach is that although the 
vibrations of the whole system are described by a con- 
tinuum of modes, the nanoparticle doesn't move as much 
for matrix phonons of certain wavelengths. 

The inelastically scattered light is modulated by the 
frequency-dependent amplitude of the vibration within 
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the nanoparticle. The mechanism for the photon scat- 
tering process depends on the electronic details of the 
nanoparticle: In semiconductors, bound excitons; In met- 
als, resonant plasmons; In insulators, strain-optic cou- 
pling. In any case, a preliminary step is to study the am- 
plitude of nanoparticle vibrations as a function of phonon 
frequency. That is the goal of our present investigation. 

The first step is to solve the equation of motion for a 
given value of R m . Extending Takagahara's workSS to 
core-shell systems, the normal modes have displacement 
fields uimnq which are normalized with respect to Eq. Q . 
Orthogonality of the modes was also checked numerically. 

For the results reported here, the outer surface of 
the matrix has zero traction force boundary conditions. 
However, we have verified that use of zero displacement 
boundary conditions has no noticeable consequence on 
the figures presented in this paper. Eigenfrequencies are 
affected, but because we are interested in the macroscopic 
limit (R m /R p — > oo) this is not important. 

To connect CSM with the FSM, BSM and CFM mod- 
els, we compare the discrete frequency sets to the position 
and width of peaks of the mean square displacement in- 
side the nanoparticle (Eq. where ui mnq is obtained 
from CSM) as a function of frequency. 



{Ulmnq 2 )p = / 1 1 U lmnq (R) || 2 d 3 R (4) 

3^Rp Jr<R p 

Using the selection rules predicted by Duval)2i one can 
get a rough idea of what the reduced Raman spectrum 
might look like: only (SPH,/=0) and (SPH,Z=2) modes 
are Raman active for a spherical isotropic nanoparticle 
under non-resonant excitation. 

The mean square displacement (u 2 ) p that we plot in 
our figures is only one of many possible measures of the 
internal motion of the nanoparticle. If we are inter- 
ested in low-frequency inelastic light scattering, then ul- 
timately we want to calculate the scattered photon spec- 
trum. Our model provides a key ingredient for such a 
calculation in the ui mnq (R). Just as in plane wave vi- 
brational modes, second quantization can be applied to 
these modes to obtain phonon creation and annihilation 
operators^ which would be required in any model of pho- 
ton scattering. 

The (u 2 )p variation with the mode frequency is not 
affected once R m /R p exceeds several times unity. In- 
creasing this ratio results in more vibrational modes for 
a given frequency range. This increase is compensated 
by the lower amplitude associated with each mode. In 
this paper, we use R m /R p = 100 when possible to have 
continuous looking plots. Calculations with bigger ra- 
tios were performed but don't change the results. Note 
that R m / R„ = 100 means that the nanoparticle occupies 
1/1000000^ of the total volume. 
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TABLE II: Wavenumber (v) and HWHM damping (Au) in 
cm - for free (FSM), BaO-P2 0s embedded Ag nanoparti- 
cles (CFM) and rigid nanoparticle (RCFM) for modes with 
v <30 cm -1 and R p = 4.9 nm. 



V. DISCUSSION 
A. Ag in BaO-P2 5 glass 

Low-frequency inelastic light scattering^ and fem- 
tosecond pump-probe experiments^ on nearly spherical 
Ag nanoparticles grown in a BaO-P205 (high density 
glass) matrix show very well resolved features, making 
it a perfect system to start with. 

Elastic parameters for silver at 300KM are direction- 
ally averaged^ giving sound speeds 1740 m/s and 3750 
m/s. The actual transverse speed of sound in Ag varies 
from 1195 m/s to 2080 m/s depending on the direction. 
Longitudinal speed of sound varies from 3410 to 3935 
m/s. Rms averaged speeds (1760 m/s, 3750 m/s) are 
comparable to the average speeds. FSM vibrational fre- 
quencies calculated using directionally averaged speeds 
are in good agreement with frequencies calculated with 
a molecular dynamic approach which fully takes into 
account elastic anisotropy^ Mass density and sound 
speeds for the BaO-P20s matrix are taken from Voisin 
et al& 
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FIG. 1: Mean square displacement within the nanoparticle 
interior of vibrational modes for Ag 7? p =4.9 nm in BaO-P20s 
with R m /R p = 100. SPH 2=0,1 and 2 and TOR 1=1 and 2 
from bottom to top. Arrows indicate positions and FWHM 
obtained with the CFM. 



to zero displacement doesn't change the appearance of 
these plots. 

Points in the (SPH,/^0) plots clearly tend to clus- 
ter on two curves, with a few points lying in between. 
The (SPH,^0) plots consist of a longitudinal-like curve 
with one broad peak and a transverse-like curve with five 
peaks. CFM calculations also reveal this feature because 
some modes are broader. 

In agreement with Voisin et al.^ Fig. ^ shows that Ag 
in BaO-P205 has significant damping, but the nanopar- 
ticle's damped pseudomodes are clearly visible. This 
probably explains why well resolved Raman peaks can 
be observed i22 

With the exception of "matrix modes" discussed in 
the following section, there is good agreement between 
the maxima of (u 2 ) p from CSM and the wavenumbers 
and HWHM as predicted using CFM. Where CFM pre- 
dicts modes where no (u 2 ) p peaks appear, these also cor- 
respond to RCFM modes as discussed below. This and 
subsequent figures demonstrate the utility of CFM to pre- 
dict the dominant frequencies at which the nanoparticle 
oscillates as well as their width. 



Table HU shows FSM, CFM and RCFM wavenumbers 
and HWHM damping. As already reported for silicon 
nanoparticles inside silicap!^ there is good agreement be- 
tween CFM and FSM. The main differences of CFM are: 
(1) nonzero damping and (2) "extra modes" as described 
in the following section. RCFM modes match some of 
the CFM modes not predicted by FSM. 

To present CSM results, the measure of nanoparti- 
cle motion in Figs. ^ to [5] is (u 2 ) p , the mean squared 
displacement within the nanoparticle of the normalized 
vibrational mode displacement field as calculated using 
Eq. J3J. In each case, (u 2 ) p is plotted individually for 
each mode versus the mode's wavenumber in cm -1 . To 
make it easier to compare different materials and differ- 
ent sizes, the different (u 2 ) p plots cover approximately 
the first two (SPH,Z=0) pseudomodes. 

Figure ^ shows individual plots for the separate cases 
of q and I. Wavenumbers and HWHM of pseudomodes 
using CFM as listed in Tab. ITU are also indicated on the 
figure as arrows. 

It is apparent that the interior of the nanoparticle is 
in motion for every CSM mode, unlike FSM wherein the 
nanoparticle would only vibrate at a discrete set of fre- 
quencies. 

Relatively smooth dependence of (u 2 ) p on frequency is 
seen for (TOR,Z=l), (TOR,Z=2) and (SPH,/=0) modes. 
The graph is much more complicated for the (SPH,/=1) 
and (SPH,Z=2) cases. In these latter spheroidal modes, 
the motion is a mixture of longitudinal and transverse 
waves travelling at different speeds. Normal modes are 
forced to mix these two polarizations in order to satisfy 
boundary conditions at R — R p and R — R m . Changing 
the boundary condition at R = R m from zero traction 



B. Extra Modes 

The presence of the matrix leads to the appearance 
of some qualitatively new modes not given by FSM. We 
classify these as: (1) semi-rigid librational modes, which 
appear at the low frequency end of the (TOR,/=l) plots 
(near 5 cm" 1 in Fig. [TJ. Both CFM and CSM con- 
firm their presence. The zero frequency limit of FSM for 
(TOR,/=l) corresponds to rigid rotation. We therefore 
identify low frequency (TOR,/=l) modes as nearly rigid 
rotational oscillation. (2) semi-rigid rattling modes, also 
at the low frequency end of the (SPH,/=1) plots (near 
3 cm -1 in Fig.UJ, also confirmed by both CFM and CSM. 
The zero frequency limit of FSM for (SPH,Z=1) corre- 
sponds to rigid rectilinear motion. We therefore identify 
low frequency (SPH,Z=1) modes as nearly rigid transla- 
tional oscillation. (3) "matrix modes" predicted by CFM 
but not apparent in the plots of (u 2 ) p from CSM. These 
modes are well approximated by RCFM, which confirms 
their nature as motions of the matrix which do not signif- 
icantly shake the nanoparticle. VCFM is a good approx- 
imation for the case of a heavier and more rigid matrix. 

CFM roots were obtained numerically in this work. 
However, some exact calculations are possible for "ma- 
trix modes". As an example, for (SPH,/=0) modes, 
there is only one VCFM eigenfrequency given by £ — 

27 (y/l - 7 2 + ijj (see Eq. (14") in Dubrovskiy et alJ£) 
and only one RCFM eigenfrequency given by £ = i where 
7 = vip/v™ and £ = ujRjv™. To our knowledge, this is 
the first report of Re(w)=0 modes in such systems. 
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C. Pseudomode Damping 

In the absence of the matrix, the nanoparticle will vi- 
brate at frequencies wfJl^. Coupling of the nanoparticle 
to the matrix clamps these pseudomodes. 

One way to get a qualitative feeling for this problem 
is to compare acoustic impedances. For plane sound 
waves in a fluid, acoustic impedance is density times 
speed of sound. Reflection at normal incidence from a 
planar interface is related to the difference of acoustic 
impedance between the two fluids. This concept also ap- 
plies in general terms to acoustic waves in a solid (see Ap- 
pendix Q. When acoustic impedances for the nanopar- 
ticle Z p = p p v p L and the matrix Z m = p m v™ match, 
reflection of longitudinal sound waves is negligible. 

As the acoustic impedance mismatch increases, acous- 
tic wave reflection increases leading to more pronounced 
nanoparticle modes. This is clearly seen in Fig. where 
varying only the density of the matrix changes the loca- 
tions of frequencies at which (u 2 ) p peaks. For low matrix 
mass density (Z p > Z m ), (u 2 ) p peaks close to the FSM 
frequencies. For high matrix mass density (Z p < Z m ), 
peaks appear close to the BSM frequencies. For matching 
impedances (p m rj 2.24pBaO-P 2 o 5 ), oscillations are weak. 
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FIG. 2: Mean square displacement within the nanoparticle 
interior of (SPH,Z=0) CSM modes for a 4.9 nm radius silver 
nanoparticle and R m /R p = 100. Matrix sound speeds are for 
BaO-P20s. Matrix mass density is scaled from BaO-P20s by 
factors of 0.5, 1, 2.24, 4 and 6 from bottom to top. Arrows 
indicate the vibrational wavenumbers for the same silver par- 
ticle with rigid (BSM, top) and free (FSM, bottom) boundary 
conditions. 



D. Na in BaO-P 2 5 glass 

This hypothetical case illustrates a soft nanoparticle 
in a hard matrix. Sodium is a very light and soft metal. 
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TABLE III: Wavenumber (u) and damping (Av) in cm 1 of 
modes of a BaO-P2 0s matrix surrounding a spherical cavity 
(VCFM) with radius R p = 4.9 nm. 



BaO-P205 is a heavy and rigid kind of glass. 

Elastic parameters for Na at 300 K are from P. Ho 
et alw^ and directionally averaged to get sound speeds 
1610 m/s and 3300 m/s. In this case Z p /Z m « 0.182. 

Figure El shows the plots of {u 2 ) p from CSM. CFM 
mode wavenumbers and HWHM are shown as arrows. 
Peaks of (u 2 ) p are always narrow and close to BSM fre- 
quencies (not shown). CFM predicts additional much 
broader peaks which do not appear in (u 2 ) p , but these 
are close to VCFM frequencies as given in Tab. IIIII Un- 
like a heavy nanoparticle in a light matrix, no additional 
low frequency mode corresponding to semi-rigid libra- 
tional (TOR,l) or rattling (SPH,1) mode appears in the 
(u 2 ) p plots. 



- — j — i 


4- 


4- 


T 


i ..^J , V... 









—I -t 




-1- 


-1- 


' ^ 


...J V... 


.J V, 







4- 


4- -1 


L „,..,, lu mmffS^L^J^^^^. 


^ 





= -HI + - 


4- 




4- 






■ 



5 10 15 20 25 30 

v (cm" ) 

FIG. 3: Mean square displacement within the nanoparticle 
interior of CSM vibrational modes for Na in BaO-P2 0s with 
R m /Rp = 100 and R p = 4.9 nm. SPH 1=0,1 and 2 and TOR 
1=1 and 2 from bottom to top. Arrows indicate positions and 
FWHM obtained from CFM. 
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TABLE IV: Wavenumber (u) and FWHM (Av) in cm" 1 of 
a matrix surrounding an infinitely rigid and heavy sphere 
(RCFM). Sphere radius is 3.4 nm. 



E. Si in Si0 2 

This example is chosen to show what happens when the 
elastic properties of the nanoparticle and the matrix are 
very similar. Low frequency Raman experiments on such 
systemsi^Si2& reveal broad peaks whose Raman position 
scales inversely with nanoparticle size. 

Elastic parameters and mass density for silicon and 
silica are taken from our previous work^ R p = 3.4 nm 
comes from the same paper and we chose to use the same 
value for the following systems. In contrast to Ag in BaO- 
P2O5 where Z p /Z m w 2.34, in the case of silicon inside 
silica, we have Z p /Z m « 1.60. As can be seen in Fig. ^ 
this results in blurred features. Silicon has a lower sound 
speed ratio (vl/vt ~ 1.67) than silver (2.15). Therefore 
fewer TOR and (SPH,Z^0) pseudomodes are seen in the 
figure. In a previous work^ we attributed the lowest 
observed low-frequency inelastic light scattering peak to 
scattering by the first (TOR,Z=l) pseudomode. At the 
same time, we pointed out that the calculated damping 
was too high. From the (TOR,Z=l) plot, we can see that 
a proper description of the scattering process is necessary 
to model the Raman spectrum. In particular, the role of 
the electronic levels need to be taken into account. At 
the same time, a better description of the silicon-silica 
interface could somewhat mechanically isolate the sili- 
con nanoparticle which would result in narrower pseu- 
domodes. Introducing a thin intermediate layer between 
the nanoparticle and the matrix is a way to achieve this. 15 
However, choosing the mass density and elastic parame- 
ters for this intermediate shell is a difficult task. 



F. CdS 

Many low-frequency Raman studies have been pub- 
lished for matrix embedded CdS^Sei-^ nanoparticles. 
Because of the different matrices used and varying 
nanoparticle composition, it is very difficult to address 
all of them here. In this paper, we will focus on CdS 
nanoparticles and show the influence of two different ma- 
trices. Low-frequency Raman spectra of CdS nanopar- 
ticles embedded in sol-gel silica have been reported by 
Othmani et alml and Saviot et <dm* Ge02 embedded CdS 
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FIG. 4: Mean square displacement within the nanoparti- 
cle interior of CSM vibrational modes for Si in Si0 2 with 
R m /R P = 100 and R p = 3.4 nm. SPH Z=0,1 and 2 and TOR 
1=1 and 2 from bottom to top. Arrows indicate positions and 
FWHM obtained from CFM. 



nanoparticles were investigated by Tanaka et alM 

Elastic parameters for CdS at 300 K are taken from 
Berlincourt et alM* and directionally-averaged as before 
to get sound speeds 1870 m/s and 4290 m/s. We use elas- 
tic constants measured at constant electric field and do 
not consider piezoelectric corrections. Mass density and 
sound speeds for amorphous Ge02 are from Antoniou et 

CSM (u 2 ) p plots and CFM results are presented in 
Fig. EJ In both cases Z p /Z m ss 1.6 which is very close to 
the previous Si in Si02 system. Ge02 having a smaller 
impedance than Si02 (even when using transverse sound 
speeds), one would expect the features of the Ge02 em- 
bedded nanoparticles to be slightly more pronounced 
than the Si02 ones because the impedance mismatch is 
bigger. It turns out to be more complex. Modes having 
a strong longitudinal character follow this trend while 
those having a strong transverse character present more 
resolved low- frequency features for Si02- This is revealed 
by both the CFM and CSM calculations. Notably, in 
the Si02 matrix, both CSM and CFM show a sharp 
(TOR,Z=l) rigid librational mode at 10.3 cm" 1 . The 
hexagonal crystal structure and lack of inversion sym- 
metry of CdS makes its electronic properties anisotropic 
even if the nanocrystal is a perfect sphere. So it is pos- 
sible that this mode could be Raman active. Acoustic 
properties are anisotropic too, even if this paper works 
around this issue by performing a directional average. 
Rigid librations may therefore have the ability to couple 
to electronic degrees of freedom and modulate inelastic 
light scattering. 

The only satisfactory way to explain the low-frequency 
Raman features is to couple the CSM approach with the 
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FIG. 5: Mean square displacement within the nanoparticle 
interior of vibrational modes for CdS in SiC>2 (top) and GeC>2 
(bottom) with R m /R P = 100 and R p = 3.4 nm. SPH 1=0,1 
and 2 and TOR 1=1 and 2 from bottom to top. Arrows indi- 
cate CFM positions and FWHM. 



of a silver nanoparticle embedded in a BaO-P20s matrix. 
The lowest frequency for a given mode (q,l) increases 
with I. The numbers of TOR and SPH modes up to 
1=60 had to be summed taking into account the 21+1 
degeneracy. In this way, the VDOS up to v = 15 cm -1 for 
Rm/Rp = 10 was obtained. Calculations for R m /R p = 3 
were also performed. Results are compared to the Debye 
model using matrix sound speeds (Eq. ©)• 

AWM _ J_ f± , 1W 

V 2tt 2 \vl + vl) 3 1 ' 

The presence of the nanoparticle doesn't appreciably 
modify the VDOS, even for small R m . 




1 10 

v (cm -1 ) 



work done by Sirenko et aim* on the interaction of the 
electronic excitations quantized in the confining potential 
of the nanoparticle with bulk- like acoustic phonons. The 
calculated Raman spectra should be further narrowed by 
taking into account the variation of the vibration ampli- 
tude inside the nanoparticle as a function of the vibration 
frequency. 

G. Vibrational density of states 



FIG. 6: Cumulative number of phonon modes per unit vol- 
ume as a function of wavenumber for a silver nanoparti- 
cle inside BaO-P205 with R m /R p = 3 (plus symbols) and 
Rm/Rp = 10 (square symbols). The corresponding matrix 
only points are connected with lines and are very close to the 
previous symbols. The straight line corresponds to the Debye 
model (Eq. ©). Position and width of the nanoparticle CFM 
pseudomodes are indicated by arrows: SPH £=0,1,2,3,4,5 and 
TOR £=1,2,3,4,5 from bottom to top respectively. "Matrix 
modes" are not shown. 



Free 3 ^ and weakly coupled 33 nanoparticles are inter- 
esting objects because their vibrational density of states 
(VDOS) below the lowest FSM mode vanish. However, 
in our case both FSM and CFM have discrete mode spec- 
tra and do not reproduce the correct phonon density of 
states. Using CSM we calculated the VDOS. Note that 
this only involves finding the frequencies of the core-shell 
system. Amplitude normalization is not needed here. 

where V is the volume of the core-shell system. 

The cumulative number of modes per unit volume up 
to wavenumber v (Eq. JSJ) is plotted in Fig.[S]for the case 



VI. CONCLUSION 

Within a continuum elastic approximation, the phonon 
spectrum of an isotropic spherical nanoparticle embedded 
in a macroscopic matrix can be calculated without ap- 
proximation. This facilitates computations that require 
information about the phonon modes, and resolves the 
apparent contradiction involved in using free sphere vi- 
brational modes to describe the mechanical oscillations 
of embedded nanoparticles. 

Even though by necessity an object embedded in a 
macroscopic matrix must vibrate at every frequency, the 
vibration of a nanoparticle in a light and soft matrix is 
clearly dominated by FSM modes. This explains the suc- 
cess of FSM to reproduce the results of Raman scattering 
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experiments on metal and semiconductor nanoparticles 
embedded in glasses. 

CFM is demonstrated to accurately predict both the 
frequency and HWHM of pseudomodes of a matrix- 
embedded nanoparticle. New low frequency features cor- 
responding to semi-rigid librational and rattling motions 
as a result of the restorative force of the matrix are clearly 
confirmed using both CFM and CSM. This points out the 
most serious limitation of FSM for embedded nanopar- 
ticles. Even so, FSM-predicted frequencies of pseudo- 
modes turn out to be surprisingly accurate. It is the 
acoustic impedance mismatch between nanoparticle and 
matrix due both to density and speed of sound differ- 
ence that allows dampening to be sufficiently slight that 
nanoparticles can " look free" in some light scattering ex- 
periments. 

While the new (SPH,Z=1) and (TOR,Z=l) modes are 
not Raman active for an isotropic perfectly spherical 
nanoparticle, various imperfections can potentially al- 
low these motions to modulate the polarizability of the 
nanoparticle and contribute to Raman and Brillouin 
spectra. However, in situations where elaboration meth- 
ods allow relatively spherical nanoparticle growth, these 
new modes remain "invisible" in Raman and Brillouin ex- 
periments. " Matrix modes" are also present, but because 
they are mainly restricted to the matrix, they can't be 
observed in experiments which are very often made under 
resonant conditions on the nanoparticle electronic tran- 
sitions. Thus, the new qualitative features of the motion 
of embedded nanoparticles are not readily observable in 
experiments. 

Size variation within a sample of nanoparticles causes 
a distribution of nanoparticle pseudomode frequencies. 
Sample characterization methods such as TEM have the 
ability to measure the extent of this size variation, but 
elaboration methods are normally limited in their abil- 
ity to restrict the range of nanoparticle radius within 
a sample. This important effect "masks" whatever line 
broadening of the pseudomodes results from matrix em- 
bedding. 

Taken together, every feature of these systems con- 
spires to obscure experimental indications that the 
nanoparticles really are embedded in solid material as 
opposed to floating in a vacuum. Consequently the orig- 
inal 1882 Lamb FSM remains a powerful calculational 
tool for the description of nanoparticle vibrations. 

Displacement fields using CFM are complex valued and 
also blow up exponentially with R. Also, the frequencies 
are complex valued. So these cannot be used as classi- 
cal degrees of freedom to be quantized. CSM provides 
a rigorous foundation for more detailed theoretical stud- 
ies since it is clear how it can be quantized. Work is in 
progress to apply this approach together with the proper 
description of the nanoparticle confined electronic levels 
in order to simulate the Raman spectrum of some sys- 
tems. 



APPENDIX A: MODE ORTHOGONALITY 

The kinetic energy and potential energy of the system 
(assuming only mode A is excited) are: 

T A (t) = j ~p(R)\\ii A (R,t)\\ 2 d 3 R (Al) 



Vk(*) = / ^CijkiUAi,jUAk,id 3 R (A2) 

By conservation of energy E A = T^t) + V A {t) cannot 
depend on time. Suppose now that two modes A and B 
(lua u>b) are simultaneously excited, so that 

u(R,t) = sin(ujAt)uA(R) + sin(u>Bt)uB{R) (A3) 

TA+B(t) is the kinetic energy due to the superposition 
of the two modes. Therefore, 

T A+B (t) = T A (t)+T B {t) (A4) 

+ U A^ 'BCOs(ui At)cos(uj gt) J pUA ' Ub^FR 

Va+b($) is the elastic potential energy due to the su- 
perposition of the two modes. So, 

V A+B (t) = V A (t) + V B (t) (A5) 
-I- sin(uJAt)sin(uiBt) J c^kiUAi^UBk.id 3 R 

With both modes excited, the energy is Ea+b- We 
next evaluate Ea+b — Ea — Eb- By energy conserva- 
tion (since there are no body or traction forces or energy 
dissipation mechanisms) E A , Eb and Ea+b must all be 
constant in time. Consequently, 

J p(R)u A ■ u B d 3 R = (A6) 



J c i:j kiUAi,jUBk,id 3 R = (A7) 

Equation (|A6|) was derived as equation (3.72) on page 
39 of Trallero-Giner et alM 

Each of Eq. 1A6|I and Eq. I|A7|I provides an orthog- 
onality on nondegenerate normal modes of vibrations. 
For convenience, we choose to orthonormalize the nor- 
mal modes as in Eq. (JIJ. Within a degenerate subspace, 
Gram-Schmidt orthogonalization can be used with this 
inner product to generate a complete orthonormal basis 
of normal modes. 
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APPENDIX B: ENERGY EQUIPARTITION 

Nanoparticle vibrations are very far from thermal 
equilibrium when a short pulse of laser light heats the 
nanoparticle to thousands of degrees. Conversely, many 
experiments are done with low incident laser beam in- 
tensities so that temperature is constant throughout the 
nanoparticle and surrounding matrix. 

For a linear classical mechanical system in thermody- 
namic equilibrium, each degree of freedom has average 
energy fceT. Here we consider the approximate applica- 
bility of this equipartition theorem to vibrations of the 
embedded nanoparticle. 

The nanoparticle is not an isolated system because of 
its mechanical connection to the matrix, but if the den- 
sity or speeds of sound are significantly different (between 
nanoparticle and matrix) then there may be in effect a 
poor impedance match at the interface so that vibrational 
energy cannot readily enter or exit the nanoparticle. In 
this case there is hope that the nanoparticle can be ap- 
proximated as a nearly isolated mechanical system only 
weakly coupled to the matrix. 

Removing the matrix, the normal vibrational "pseu- 
domodes" of the nanoparticle can be indexed by Imnq. 
The term " pseudomodes" emphasizes that they are not 
normal modes because of the non-isolation of the sys- 
tem. Three indices I, m and q are in common with the 
corresponding indices of modes of the core-shell system 
as a whole, n is the index for the nanoparticle pseudo- 
modes and n is the index for normal modes of the entire 
core-shell system. 

Assuming weak coupling, each pseudomode Imhq will 
have average kinetic energy ifcgT. The coupling to the 
matrix leads to thermal fluctuations both of amplitude 
and phase of the motion within the nanoparticle. As a 



result the motion of a single pseudomode has a widened 
band of frequencies. 

We add up kinetic energies within the nanoparticle for 
all core-shell modes whose u>i mnq is close to u>i m hq of the 
pseudomode. This gives a time average kinetic energy 
associated with the pseudomode: 

{Tlmhq)t = ^hnnq( X f7ji7iq+yhnnq)Mp( U hnnq)p (Bl) 

n 

Since each core-shell mode has average energy fc^T, 

tfmnq+yfmnq ec l uals 2k B T 7 ' ( M P+™uf mnq ) ■ Thus, 

5>U>* « (B2) 

n ™ 

where the summation over n includes only such core-shell 
modes for which w; mrig is close to a particular pseudo- 
mode frequency cuimnq ■ This is well satisfied for the cases 
we have checked. 

Since normal modes are normalized with respect to 
Eq. Q and since the integrand of (0} is non-negative, 

(uLnq) P < M p+m /M p . 

APPENDIX C: ACOUSTIC IMPEDANCE 

Acoustic impedance for sound waves can be general- 
ized to the situation of a spherical wave crossing a spher- 
ical surface at normal incidence. In this case the acous- 
tic impedance is complex valued. In addition to density 
and speed of sound it also depends on frequency and 
the radius of the interface. When the interface radius is 
large compared to the wavelength the acoustic impedance 
again approaches density times speed of sound. See equa- 
tion (5.11.12) on page 128 of L. E. Kinsler et alM 
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